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Abstract 



The a-oj model of nuclei is studied at leading order in the expansion 
thereby introducing the self consistent Hartree approximation, the Dirac sea 
corrections and the one fermion loop meson self energies in a unified way. For 
simplicity, the Dirac sea is further treated within a semiclassical expansion 
to all orders. The well-known Landau pole vacuum instability appearing in 
this kind of theories is removed by means of a scheme recently proposed 
in this context. The effect of such removal on the low momentum effective 
parameters of the model, relevant to describe nuclear matter and finite nuclei, 
is analyzed. The one fermion loop meson self energies are found to have a 
sizeable contribution to these parameters. However, such contribution turns 
out to come mostly from the Landau poles and is thus spurious. We conclude 
that the fermionic loop can only be introduced consistently in the a-u} nuclear 
model if the Landau pole problem is dealt with properly. 
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I. INTRODUCTION 



The relativistic approach to nuclear physics has attracted much attention. From a the- 
oretical point of view, it allows to implement, in principle, the important requirements of 
relativity, unitarity, causality and renormalizability From the phenomenological side, it 
has also been successful in reproducing a large body of experimental data . In the con- 
text of finite nuclei a large amount of work has been done at the Hartree level but considering 
only the positive energy single particle nucleon states. The Dirac sea has also been studied 
since it is required to preserve the unitarity of the theory. Actually, Dirac sea corrections 
have been found to be non negligible using a semiclassical expansion which, if computed to 
fourth order, seems to be quickly convergent [^. Therefore, it would appear that the overall 
theoretical and phenomenological picture suggested by the relativistic approach is rather 
reliable. 

However, it has been known since ten years that such a description is internally incon- 
sistent. The vacuum of the theory is unstable due to the existence of tachyonic poles in the 
meson propagators at high Euclidean momenta 0. Alternatively, a translationally invariant 
mean field vacuum does not correspond to a minimum; the Dirac sea vacuum energy can 
be lowered by allowing small size mean field solutions . Being a short distance instability 
it does not show up for finite nuclei at the one fermion loop level and within a semiclassi- 
cal expansion (which is an asymptotic large size expansion). For the same reason, it does 
not appear either in the study of nuclear matter if translational invariance is imposed as a 
constraint. However, the instability sets in either in an exact mean field valence plus sea 
(i.e., one fermion loop) calculation for finite nuclei or in the determination of the correla- 
tion energy for nuclear matter (i.e., one fermion loop plus a boson loop). Unlike quantum 
electrodynamics, where the instability takes place far beyond its domain of applicability, in 
quantum hadrodynamics it occurs at the length scale of 0.2 fm that is comparable to the 
nucleon size and mass. Therefore, the existence of the instability contradicts the original 
motivation that lead to the introduction of the field theoretical model itself. In such a 
situation several possibilities arise. Firstly, one may argue that the model is defined only 
as an effective theory, subjected to inherent limitations regarding the Dirac sea. Namely, 
the sea may at best be handled semiclassically, hence reducing the scope of applicability of 
the model. This interpretation is intellectually unsatisfactory since the semiclassical treat- 
ment would be an approximation to an inexistent mean field description. Alternatively, 
and taking into account the phenomenological success of the model, one may take more 
seriously the spirit of the original proposal |Q, namely, to use specific renormalizable La- 
grangians where the basic degrees of freedom are represented by nucleon and meson fields. 
Such a path has been explored in a series of papers [P|-pT|] inspired by the early work of 
Redmond and Bogolyubov on non asymptotically free theories [ll2| , p!3| . The key feature of 
this kind of theories is that they are only defined in a perturbative sense. According to the 
latter authors, it is possible to supplement the theory with a prescription based on an exact 
fulfillment of the Kallen-Lehmann representation of the two point Green's functions. The 
interesting aspect of this proposal is that the Landau poles are removed in such a way that 
the perturbative content of the theory remains unchanged. In particular, this guarantees 
that the perturbative renormalizability is preserved. It is, however, not clear whether this 
result can be generalized to three and higher point Green's functions in order to end up with 
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a completely well-behaved field theory. Although the prescription to eliminate the ghosts 
may seem to be ad hoc, it certainly agrees more with the original proposal and provides a 
workable calculational scheme. 

The above mentioned prescription has already been used in the context of nuclear physics. 
In ref. ||T^, it was applied to ghost removal in the o exchange in the iVA^ potential. More 
recently, it has been explored to study the correlation energy in nuclear matter in the a-io 
model [Tl[] and also in the evaluation of response functions within a local density approxi- 
mation 0. Although this model is rather simple, it embodies the essential field theoretical 
aspects of the problem while still providing a reasonable phenomenological description. We 
will use the a-io model in the present work, to estimate the binding energy of finite nuclei 
within a self-consistent mean field description, including the effects due to the Dirac sea, 
after explicit elimination of the ghosts. An exact mean field calculation, both for the valence 
and sea, does make sense in the absence of a vacuum instability but in practice it becomes 
a technically cumbersome problem. This is due to the presence of a considerable number of 
negative energy bound states in addition to the continuum states 0. Therefore, it seems 
advisable to use a simpler computational scheme to obtain a numerical estimate. This will 
allow us to see whether or not the elimination of the ghosts induces dramatic changes in 
the already satisfactory description of nuclear properties. In this work we choose to keep 
the full Hartree equations for the valence part but employ a semiclassical approximation 
for the Dirac sea. This is in fact the standard procedure As already mentioned, and 

discussed in previous work 0, this expansion converges rather quickly and therefore might 
be reliably used to estimate the sea energy up to possible corrections due to shell effects. 

The paper is organized as follows. In section ^ we present the a-io model of nuclei in the 
1/A^ leading approximation, the semiclassical treatment of the Dirac sea, the renormalization 
prescriptions and the different parameter fixing schemes that we will consider. In section |T| 
we discuss the vacuum instability problem of the model and Redmond's proposal. We also 
study the implications of the ghost subtraction on the low momentum effective parameters. 
In section |V| we present our numerical results for the parameters as well as binding energies 
and mean quadratic charge radii of some closed-shell nuclei. Our conclusions are presented 
in section |V[ Explicit expressions for the zero momentum renormalized meson self energies 
and related formulas are given in the appendix. 



II. a-uj MODEL OF NUCLEI 

In this section we revise the o-lo model description of finite nuclei disregarding throughout 
the instability problem; this will be considered in the next section. The Dirac sea corrections 
are included at the semiclassical level and renormalization issues as well as the various ways 
of fixing the parameters of the model are also discussed here. 

A. Field theoretical model 

Our starting point is the Lagrangian density of the g-uj model [||,|]-|^ given by 

L{x) = ^(x) [^^{td^ - g.V^ix)) - (M - gs<p{x))] ^{x) + ^ {d^<f){x)dy{x) - <p\x)) 
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-- + - mlV^{x)V^{x) + 5C{x) . (1) 

is the isospinor nucleoli field, 0(x) the scalar field, V^(a;) the cj-meson field and 
Ffj^u = dfjYi, — dt,V^. In the former expression the necessary counterterms required by 
renormalization are accounted for by the extra Lagrangian term 6C{x) (including meson 
self-couplings). 

Including Dirac sea corrections requires to take care of renormalization issues. The best 
way of doing this in the present context is to use an effective action formalism. Further we 
have to specify the approximation scheme. The effective action will be computed at lowest 
order in the 1/N expansion, being the number of nucleon species (with Qs and of order 
l/-\/iV), that is, up to one fermion loop and tree level for bosons |TT|. This corresponds to 
the Hartree approximation for fermions including the Dirac sea p^ . 

In principle, the full effective action would have to be computed by introducing bosonic 
and fermionic sources. However, since we will consider only stationary situations, we do not 
need to introduce fermionic sources. Instead, we will proceed as usual by integrating out 
exactly the fermionic degrees of freedom. This gives directly the bosonic effective action at 
leading order in the expansion: 

T[^,V] = Tb[^,V] + Tf[(^,V], (2) 

where 

rB[0, y]=J{l i9,<l>dy - ml 02) _ 1 F^^F^'^ + i mlV.V^^ d'x , (3) 

and 

r^[0, V] = -I log Det [7^(29'^ - g.y^') - (M - g,<P)] + j 6C{x)d'x . (4) 

The fermionic determinant can be computed perturbatively, by adding up the one-fermion 
loop amputated graphs with any number of bosonic legs, using a gradient expansion or by 
any other technique. The ultraviolet divergences are to be canceled with the counterterms 
by using any renormalization scheme; all of them give the same result after fitting to physical 
observables. 

The effective action so obtained is uniquely defined and completely finite. However, 
there still remains the freedom to choose different variables to express it. We will work 
with fields renormalized at zero momentum. That is, the bosonic fields (f){x) and V^{x) are 
normalized so that their kinetic energy term is the canonical one. This is the choice shown 
above in Tb[4>, V]. Other usual choice is the on-shell one, namely, to rescale the fields so that 
the residue of the propagator at the meson pole is unity. Note that the Lagrangian mass 
parameters and do not correspond to the physical masses (which will be denoted 
m„ and m^^ in what follows) since the latter are defined as the position of the poles in 
the corresponding propagators. The difference comes from the fermion loop self energy in 
Tp[(j), V] that contains terms quadratic in the boson fields with higher order gradients. 

Let us turn now to the fermionic contribution, Tplnj), V]. We will consider nuclear ground 
states of spherical nuclei, therefore the space-like components of the a;-meson field vanish 
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0] and the remaining fields, and Vq{x) are stationary. As it is well-known, for sta- 

tionary fields the fermionic energy, i.e., minus the action rj7[0, V^] per unit time, can be 
formally written as the sum of single particle energies of the fermion moving in the bosonic 



background |T5 



Ef[<P,Vo] = Y.E^, (5) 

n 

and 

[-ia ■ V + gvVoix) + f3{M - ^^^(a:))] i)n{x) = En i)n{x) . (6) 

Note that what we have called the fermionic energy contains not only the fermionic kinetic 
energy, but also the potential energy coming from the interaction with the bosons. 

The orbitals, and thus the fermionic energy, can be divided into valence and sea, i.e., 
positive and negative energy orbitals. In realistic cases there is a gap in the spectrum which 
makes such a separation a natural one. The valence energy is therefore given by 

E7\ct^,V] = Y.ET'. (7) 

n 

On the other hand, the sea energy is ultraviolet divergent and requires the renormalization 
mentioned above The (at zero momentum) renormalized sea energy is known in a 
gradient or semiclassical expansion up to fourth order and is given by 




Ef'^= ^ , [d^xl - ll$-^(V$)^-22<l>-^(V\/)2(V$)2 + 44<l>-^((Vi$)(Viy))^ 
57607r"' J [ 

-44<l>-3((Vi$)(Vi\^))(VV) -8$-^(Vl^)^ + 22<l>-3(v2$)(V$)2 

+14 (^-%VV)^{\/^^) - 18 $-2(v2$)2 + 24 <l>-2(vV)2| . (8) 

Here, V = gvVo, $ = M — g^cj) and 7 is the spin and isospin degeneracy of the nucleon, i.e., 
2N if there are N nucleon species (in the real world N = 2). The sea energy is obtained by 
adding up the terms above. The fourth and higher order terms are ultraviolet finite as follows 
from dimensional counting. The first two terms, being renormalized at zero momentum, do 
not contain operators with dimension four or less, such as 0^, 0^, or (VVy, since they 
are already accounted for in the bosonic term rB[0, \^]. Note that the theory has been 
renormalized so that there are no three- or four-point bosonic interactions in the effective 
action at zero momentum P]. 

By definition, the true value of the classical fields (i.e., the value in the absence of external 
sources) is to be found by minimization of the effective action or, in the stationary case, of 
the energy 

E[4>, V] = Eb[(P, V] + Ef[(P, V] + V] . (9) 
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Such minimization yields the equations of motion for the bosonic fields: 



{V'-ml)ct>{x) = -g,[pT\x)+p7\x) 

(V^ - ml)V,{x) = -g, [p^-'(a;) + p''-\x)\ . (10) 
Here, Ps{x) = (\E'(a;)^'(a;)) is the scalar density and p{x) = (\E'"''(x)\l'(x)) the baryonic one: 

1 r 77'™' (sea) 
val (sea)/ N _ ^ '^^F 

9s 0(f){x) 

1 X 7?™' (sea) 

p-^(-^)(x) = +1^9^. (11) 



The set of bosonic and fermionic equations, eqs. (|T0D and (y) respectively, are to be solved 
self-consistently. Let us remark that treating the fermionic sea energy using a gradient or 
semiclassical expansion is a further approximation on top of the mean field approximation 
since it neglects possible shell effects in the Dirac sea. However, a direct solution of the 
mean field equations including renormalization of the sum of single-particle energies would 
not give a physically acceptable solution due to the presence of Landau ghosts. They will 
be considered in the next section. 

At this point it is appropriate to make some comments on renormalization. As we have 
said, one can choose different normalizations for the mesonic fields and there are also several 
sets of mesonic masses, namely, on-shell and at zero momentum. If one were to write the 
mesonic equations of motion directly, by similarity with a classical treatment, there would 
be an ambiguity as to which set should be used. The effective action treatment makes it 
clear that the mesonic field and masses are those at zero momentum. On the other hand, 
since we have not included bosonic loops, the fermionic operators in the Lagrangian are not 
renormalized and there are no proper vertex corrections. Thus the nucleon mass M, the 
nuclear densities (^'^') and the combinations gs4'{x) and gvVfj,{x) are fixed unambiguously in 
the renormalized theory. The fermionic energy V], the potentials $(x) and V{x) and 

the nucleon single particle orbitals are all free from renormalization ambiguities at leading 
order in 



B. Fixing of the parameters 



The a-ui and related theories are effective models of nuclear interaction, and hence their 
parameters are to be fixed to experimental observables within the considered approximation. 
Several procedures to perform the fixing can be found in the literature the more 

sophisticated versions try to adjust, by minimizing the appropriate function, as many 
experimental values as possible through the whole nuclear table These methods are 
useful when the theory implements enough physical elements to provide a good description 
of atomic nuclei. The particular model we are dealing with can reproduce the main features 
of nuclear force, such as saturation and the correct magic numbers; however it lacks many of 
the important ingredients of nuclear interaction, namely Coulomb interaction and p and vr 
mesons. Therefore, we will use the simple fixing scheme proposed in ref. for this model. 
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Initially there are five free parameters: the nucleon mass (M), two boson Lagrangian 
masses (m^ and m^) and the corresponding coupling constants {qs and Qy). The five physical 
observables to be reproduced are taken to be the physical nucleon mass, the physical u- 
meson mass m^, the saturation properties of nuclear matter (binding energy per nucleon 
B/A and Fermi momentum kp) and the mean quadratic charge radius of ^''Ca. In our 
approximation, the equation of state of nuclear matter at zero temperature, and hence its 
saturation properties, depends only on the nucleon mass and on ms,v and gs,v through the 
combinations [0] 

= 9s T 1 = 9v T • (12) 

ITLg ITLy 

At this point, there still remain two parameters to be fixed, e.g., m^, and Qs- Now we 
implement the physical a;-meson mass constraint. From the expression of the uj propagator 
at the leading approximation, we can obtain the value of the physical oj pole as a 
function of the Lagrangian parameters M, and or more conveniently as a function 
of M, Cy and (see appendix). Identifying the uj pole and the physical u) mass, and 
given that M and C„ have already been fixed, we obtain the value of m^. Finally, the value 
of Qs is adjusted to fit the mean quadratic charge radius of "^^Ca. We will refer to this 
fixing procedure as the u-shell scheme: the name stresses the correct association between 
the pole of the tu-meson propagator and the physical uj mass. The above fixing procedure 
gives different values of and Qs depending on the order at which the Dirac sea energy is 
included in the semiclassical expansion (see section |rV|). 

Throughout the literature the standard fixing procedure when the Dirac sea is included 
has been to give to the Lagrangian mass m„ the value of the physical uj mass (see, 
however, refs. |10,^). Of course, this yields a wrong value for the position of the u;- meson 



propagator pole, which is underestimated. We will refer to this procedure as the naive 
scheme. Note that when the Dirac sea is not included at all, the right viewpoint is to 
consider the theory at tree level, and the cu-shell and the naive schemes coincide. 

III. LANDAU INSTABILITY SUBTRACTION 

As already mentioned, the a-uj model, and more generally any Lagrangian which couples 
bosons with fermions by means of a Yukawa-like coupling, exhibits a vacuum instability . 
This instability prevents the actual calculation of physical quantities beyond the mean field 
valence approximation in a systematic way. Recently, however, a proposal by Redmond 



that explicitly eliminates the Landau ghost has been implemented to describe relativistic 
nuclear matter in a series of papers P|-pT]]. The main features of such method are contained 
already in the original papers and many details have also been discussed. For the sake of 
clarity, we outline here the method as applies to the calculation of Dirac sea effects for 
closed- shell finite nuclei. 



A. Landau instability 

Since the Landau instability shows up already at zero nuclear density, we will begin 
by considering the vacuum of the a-uj theory. On a very general basis, namely, Poincare 
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invariance, unitarity, causality and uniqueness of the vacuum state, one can show that the 
two point Green's function (time ordered product) for a scalar field admits the Kallen- 
Lehmann representation |T6| 




(13) 



where the full propagator in the vacuum is 



D{x' — x) 



^(O|T0(xO0(x)|O), 



(14) 



and the free propagator reads 



Dq{x' — X ; /i^) 



/ 



(15) 



The spectral density p(/i^) is defined as 



p{q') = {2'nfY.^\pr,-qmm\n)\\ 



(16) 



n 



It is non negative, Lorentz invariant and vanishes for space-like four momentum q. 

The Kallen-Lehmann representability is a necessary condition for any acceptable theory, 
yet it is violated by the a-uj model when the meson propagators are approximated by their 
leading term. It is not clear whether this failure is tied to the theory itself or it is an 
artifact of the approximation — it is well-known that approximations to the full propagator 
do not necessarily preserve the Kallen-Lehmann representability — . The former possibility 
would suppose a serious obstacle for the theory to be a reliable one. 

In the above mentioned approximation, eq. (|13|) still holds both for the a and the uj cases 
(in the latter case with obvious modification to account for the Lorentz structure) but the 
spectral density gets modified to be 



where p^'"(yU^) is a physically acceptable spectral density, satisfying the general requirements 
of a quantum field theory. On the other hand, however, the extra term spoils these general 
principles. The residue —Rg is negative, thus indicating the appearance of a Landau ghost 
state which contradicts the usual quantum mechanical probabilistic interpretation. More- 
over, the delta function is located at the space-like squared four momentum —Mq indicating 
the occurrence of a tachyonic instability. As a perturbative analysis shows, the dependence 
of Rg and Mg with the fermion-meson coupling constant g in the weak coupling regime is 
Rg ~ and Mq ~ 4M^ exp(47r^/(yf^), with M the nucleon mass. Therefore the perturba- 
tive content of p(/U^) and p^^i^fi"^) is the same, i.e., both quantities coincide order by order 
in a power series expansion of g keeping /x^ fixed. This can also be seen in the propagator 
form of the previous equation 




(17) 



D{p)=D^'^{p) 



Rg 



(18) 



p2 + ■ 
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For fixed four momentum, tlie gliost term vanislies as exp(— 47r^/(7^) wlien tlie coupling 
constant goes to zero. As noted by Redmond it is tlierefore possible to modify the 
theory by adding a suitable counterterm to the action that exactly cancels the ghost term 
in the meson propagator without changing the perturbative content of the theory. In this 
way the full meson propagator becomes D^^{p) which is physically acceptable and free from 
vacuum instability at leading order in the expansion. 

It is not presently known whether the stability problems of the original a-u theory are 
intrinsic or due to the approximation used, thus Redmond's procedure can be interpreted 
either as a fundamental change of the theory or as a modification of the approximation 
scheme. Although both interpretations use the perturbative expansion as a constraint, it is 
not possible, at the present stage, to decide between them. It should be made quite clear that 
in spite of the seemingly arbitrariness of the no-ghost prescription, the original theory itself 
was ambiguous regarding its non perturbative regime. In fact, being a non asymptotically 
free theory, it is not obvious how to define it beyond finite order perturbation theory. For 
the same reason, it is not Borel summable and hence additional prescriptions are required 
to reconstruct the Green's functions from perturbation theory to all orders. As an example, 
if the nucleon self energy is computed at leading order in a expansion, the existence 
of the Landau ghost in the meson propagator gives rise to a pole ambiguity. This is unlike 
physical time-like poles, which can be properly handled by the customary +ie rule, and 
thus an additional ad hoc prescription is needed. This ambiguity reflects in turn in the 
Borel transform of the perturbative series; the Borel transform presents a pole, known as 



renormalon in the literature [|l^. In recovering the sum of the perturbative series through 
inverse Borel transformation a prescription is then needed, and Redmond's proposal provides 
a particular suitable way of fixing such ambiguity. Nevertheless, it should be noted that even 
if Redmond's prescription turns out to be justified, there still remains the problem of how 
to extend it to the case of three- and more point Green's functions, since the corresponding 
Kallen-Lehmann representations has been less studied. 



B. Instability subtraction 

To implement Redmond's prescription in detail we start with the zero-momentum renor- 
malized propagator in terms of the proper self-energy for the scalar field (a similar construc- 
tion can be carried out for the vector field as well), 

D,{p') = ip'-ml-n,{p'))-\ (19) 

where the rris is the zero-momentum meson mass and the corresponding renormalization 
conditions are 11^(0) = n'^(O) = 0. The explicit formulas for the scalar and vector meson self 
energies are given in the appendix. Of course, DgijP') is just the inverse of the quadratic part 
of the effective action KgijP'). According to the previous section, the propagator presents a 
tachyonic pole. Since the ghost subtraction is performed at the level of the two-point Green's 
function, it is clear that the corresponding Lagrangian counterterm must involve a quadratic 
operator in the mesonic fields. The counterterm kernel AKgijP') must be such that cancels 
the ghost term in the propagator Dgip^) in eq. ([18|). The subtraction does not modify the 
position of the physical meson pole nor its residue, but it will change the zero-momentum 
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parameters and also the off-shell behavior. Both features are relevant to nuclear properties. 
This will be discussed further in the next section. 
Straightforward calculation yields 

As stated, this expression vanishes as exp{—ATT^/g^) for small Qs at fixed momentum. There- 
fore it is a genuine non perturbative counterterm. It is also non local as it depends in a 
non polynomial way on the momentum. In any case, it does not introduce new ultraviolet 
divergences at the one fermion loop level. However, it is not known whether the presence of 
this term spoils any general principle of quantum field theory. 

Proceeding in a similar way with the cu-field V^{x), the following change in the total 
original action is induced 

= i / ^,<P{-P)AK,{P')<I>{P) -\j ^V;.(-p)Airrb')K(p) , (21) 

where (j){p) and V^{p) are the Fourier transform of the scalar and vector fields in coordinate 
space, (j){x) and V^{x) respectively. Note that at tree- level for bosons, as we are considering 
throughout, this modification of the action is to be added directly to the effective action 
— in fact, this is the simplest way to derive eq. (^OD — . Therefore, in the case of static fields, 
the total mean field energy after ghost elimination reads 

E = E'f + E''^'' + Eb + AE, (22) 

where E],^\ and Eb were given in section |I| and 

AE[c^, V] = ^J d'x(P{x)AKs{V'mx) j d''xVo{x)AK'!^{V^)Vo{x) . (23) 

One can proceed by minimizing the mean field total energy as a functional of the bosonic 
and fermionic fields. This yields the usual set of Dirac equations for the fermions, eqs. 
and modifies the left-hand side of the bosonic eqs. (|TUp by adding a linear non-local term. 
This will be our starting point to study the effect of eliminating the ghosts in the description 
of finite nuclei. We note that the instability is removed at the Lagrangian level, i.e., the non- 
local counterterms are taken to be new terms of the starting Lagrangian which is then used 
to describe the vacuum, nuclear matter and finite nuclei. Therefore no new prescriptions 
are needed in addition to Redmond's to specify how the vacuum and the medium parts of 
the effective action are modified by the removal of the ghosts. 

So far, the new counterterms, although induced through the Yukawa coupling with 
fermions, have been treated as purely bosonic terms. Therefore, they do not contribute 
directly to bilinear fermionic operators such as baryonic and scalar densities. An alternative 
viewpoint would be to take them rather as fermionic terms, i.e., as a (non-local and non- 
perturbative) redefinition of the fermionic determinant. The energy functional, and thus the 
mean field equations and their solutions, coincide in the bosonic and fermionic interpreta- 
tions of the new term, but the baryonic densities and related observables would differ, since 
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they pick up a new contribution given the corresponding formulas similar to eqs. ([TTD. Am- 
biguities and redefinitions are ubiquitous in quantum field theories, due to the well-known 
ultraviolet divergences. However, in well-behaved theories the only freedom allowed in the 
definition of the fermionic determinant comes from adding counterterms which are local and 
polynomial in the fields. Since the new counterterms induced by Redmond's method are not 
of this form, we will not pursue such alternative point of view in what follows. Nevertheless, 
a more compelling argument would be needed to make a reliable choice between the two 
possibilities. 



C. Application to finite nuclei 

In this section we will take advantage of the smooth behavior of the mesonic mean fields 
in coordinate space which allows us to apply a derivative or low momentum expansion. The 
quality of the gradient expansion can be tested a posteriori by a direct computation. The 
practical implementation of this idea consists of treating the term AS by expanding each of 
the kernels AK{p'^) in a power series of the momentum squared around zero 

AKip^) = Y.AK2np'''. (24) 

n>0 

The first two terms are given explicitly by 

AK, - "^'^^ 



{Ml - rn?RG) 



The explicit expressions of the tachyonic pole parameters Mg and Rg for each meson can 
be found below. 

Numerically, we have found that the fourth and higher orders in this gradient expansion 
are negligible as compared to zeroth- and second orders. In fact, in ref. the same behavior 
was found for the correction to the Dirac sea contribution to the binding energy of a nucleus. 
As a result, even for light nuclei, Ef^^ in eq. (H) can be safely neglected. Furthermore, it has 
been shown |]TB[ that the fourth order term in the gradient expansion of the valence energy, 
if treated semiclassically, is less important than shell effects. So, it seems to be a general rule 
that, for the purpose of describing static nuclear properties, only the two lowest order terms 
of a gradient expansion need to be considered. We warn, however, that the convergence of 
the gradient or semiclassical expansion is not the same as converging to the exact mean field 
result, since there could be shell effects not accounted for by this expansion at any finite 



order. Such effects, certainly exist in the valence part [|18|. Even in a seemingly safe case 
as infinite nuclear matter, where only the zeroth order has a non vanishing contribution, 
something is left out by the gradient expansion since the exact mean field solution does not 
exist due to the Landau ghost instability (of course, the situation may change if the Landau 
pole is removed). In other words, although a gradient expansion might appear to be exact 
in the nuclear matter case, it hides the very existence of the vacuum instability. 

From the previous discussion it follows that the whole effect of the ghost subtraction 
is represented by adding a term AS* to the effective action with same form as the bosonic 
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part of the original theory, TB[(f), V] in eq. (^. This amounts to a modification of the zero- 
momentum parameters of the effective action. The new zero- momentum scalar field (i.e., 
with canonical kinetic energy), mass and coupling constant in terms of those of the original 
theory are given by 

0(x) = (l + Air2')i/2^(x), 



1/2 



= ^^-rW^ , (26) 




The new coupling constant is obtained recalling that gs4'{x) should be invariant. Similar 
formulas hold for the vector meson. With these definitions (and keeping only AKg ^lp"^) till 
second order in p"^) one findsQ 

EB[(f),V]ms,m^] = EB[4>,V;ms,my] + AE[(j),V]ms,my] , (27) 
Ef&, V; Qs, dv] = V; Qs, Qv] ■ 

The bosonic equations for the new meson fields after ghost removal are hence identical to 
those of the original theory using 

as zero-momentum masses and coupling constants respectively. In the limit of large ghost 
masses or vanishing ghost residues, the reparameterization becomes trivial, as it should be. 
Let us note that although the zero-momentum parameters of the effective action m^^^ and 
Qs^v are the relevant ones for nuclear structure properties, the parameters rus^v and g^^v are 
the (zero-momentum renormalized) Lagrangian parameters and they are also needed, since 
they are those appearing in the Feynman rules in a perturbative treatment of the model. Of 
course, both sets of parameters coincide when the ghosts are not removed or if there were 
no ghosts in the theory. 

To finish this section we give explicitly the fourth order coefficient in the gradient ex- 
pansion of AE^ taking into account the rescaling of the mesonic fields, namely, 

AK^ _ Rc^Ml^m'f ^g" Rl - Im" MIRg 

■ 5 i^yj 



l + Ai^a [M% ^ Rg) [Ml - Rg) an'^ + Rg 

where a is 160 for the scalar meson and 120 for the vector meson. As already stated, 
for typical mesonic profiles the contribution of these fourth order terms are found to be 
numerically negligible. Simple order of magnitude estimates show that squared gradients 



-"^Note that Eb^f[ ] refer to the functionals (the same at both sides of the equations) and not to 
their value as is also usual in physics literature. 



12 



are suppressed by a factor {RMg)~^, R being the nuclear radius, and therefore higher orders 
can also be neglected. That the low momentum region is the one relevant to nuclear physics 
can also be seen from the kernel ii's(p^), shown in fig. |I[ From eq. (|21|), this kernel is to be 
compared with the function (j){p) that has a width of the order of R~^. It is clear from the 
figure that at this scale all the structure of the kernel at higher momenta is irrelevant to 
AE. 



D. Fixing of the parameters after ghost subtraction 

As noted in section ||, the equation of state at zero temperature for nuclear matter 
depends only on the dimensionless quantities and C^, that now become 

C! = fs — , Cl = gl — . 30 

Fixing the saturation density and binding energy to their observed values yields, of course, 
the same numerical values for and as in the original theory. After this is done, all 
static properties of nuclear matter are determined and thus they are insensitive to the ghost 
subtraction. Therefore, at leading order in the 1/N expansion, to see any effect one should 
study either dynamical nuclear matter properties as done in ref. [|rO| or finite nuclei as we 
do here. 

It is remarkable that if all the parameters of the model were to be fixed exclusively 
by a set of nuclear structure properties, the ghost subtracted and the original theories 
would be indistinguishable regarding any other static nuclear prediction, because bosonic 
and fermionic equations of motion have the same form in both theories. They would differ 
however far from the zero four momentum region where the truncation of the ghost kernels 
AK{p^) at order is no longer justified. In practice, the predictions will change after ghost 
removal because the cj-meson mass is quite large and is one of the observables to be used in 
the fixing of the parameters. 

To fix the parameters of the theory we choose the same observables as in section |I[ Let 
us consider first the vector meson parameters and g^- We proceed as follows: 

1. We choose a trial value for (the zero- momentum coupling constant of the original 
theory). This value and the known physical values of the w-meson and nucleon masses, 
and M respectively, determines (the zero- momentum mass of the original theory), 
namely 



,2 2 



3 9 M2 3 V MV \ m2 \2M 



— 1 arcsm —— I > . (31) 



(This, as well as the formulas given below, can be deduced from those in the appendix.) 

2. g^ and provide the values of the tachyonic parameters R"q and Mq. They are 
given by 



Ml 



JDV 



2M 




(32) 
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where the quantity k„ is the real solution of the following equation (there is an imaginary 
solution which corresponds to the cu-meson pole) 

^ ' (-^-l) + ^lf|-^Vog^--^ + ?| = 0. (33) 



3. Known g^, m„, Mq and i?^, the values of and are obtained from eqs. 
They are then inserted in eqs. (|30|) to yield C^. If necessary, the initial trial value of g^ 
should be readjusted so that the value of so obtained coincides with that determined by 
the saturation properties of nuclear matter. 

The procedure to fix the parameters and gs is similar but slightly simpler since the 
physical mass of the scalar meson mg- is not used in the fit. Some trial values for m<j and gs 
are proposed. This allows to compute Mq and Rq by means of the formulas 



2M 



R'q 16 7r2 

where is the real solution of 





(34) 



^ + 4#''=--^)-iif5H'°^;^3i-2'=2-5} = "- (35) 

One can then compute rhs and gs and thus Cl and the mean quadratic charge radius of 
''^Ca. The initial values of rus and gs should be adjusted to reproduce these two quantities. 
We will refer to the set of masses and coupling constants so obtained as the no-ghost scheme 
parameters. 



IV. NUMERICAL RESULTS AND DISCUSSION 

As explained in section |I|, the parameters of the theory are fitted to five observables. 
For the latter we take the following numerical values: M = 939 MeV, = 783 MeV, 
B/A = 15.75 MeV, kp = 1.3 fm~^ and 3.82 fm for the mean quadratic charge radius of 
40Ca. 

If the Dirac sea is not included at all, the numerical values that we find for the nuclear 
matter combinations and are 

= 357.7, C^2 = 274.1 (36) 

The corresponding Lagrangian parameters are shown in table |I[ There we also show and 
that correspond to the position of the poles in the propagators after including the one- 
loop meson self energy. They are an output of the calculation and are given for illustration 
purposes. 

When the Dirac sea is included, nuclear matter properties fix the following values 

= 227.8 , = 147.5 (37) 
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Note that in nuclear matter only the zeroth order Eq^'^ is needed in the gradient expansion 
of the sea energy, since the meson fields are constant. The (zero momentum renormalized) 
Lagrangian meson masses and coupling constants Qg^y are shown in table | in various 
schemes, namely, tu-shell, no-ghost and naive schemes, previously defined. The scalar meson 
parameters differ if the Dirac sea energy is included at zeroth order or at all orders (in 
practice zeroth plus second order) in the gradient expansion. For the sake of completeness, 
both possibilities are shown in the table. The numbers in brackets in the no-ghost scheme are 
the zero-momentum parameters of the effective action, m<, and Qs^v (in the other schemes 
they coincide with the Lagrangian parameters). Again m^j and m^^ refer to the scalar and 
vector propagator-pole masses after including the one fermion loop self energy for each set 
of Lagrangian parameters. Table |T| shows the ghost masses and residues corresponding to 
the zero-momentum renormalized propagators. The no-ghost scheme parameters have been 
used. 

The binding energies per nucleon (without center of mass corrections) and mean 
quadratic charge radii (without convolution with the nucleon form factor) of several closed- 
shell nuclei are shown in tables ^ and |V| for the cu-shell and for the naive and no-ghost 
schemes (these two schemes give the same numbers), as well as for the case of not including 
the Dirac sea. The experimental data are taken from refs. p!U|-pi[]. 

From table | it follows that the zero-momentum vector meson mass in the tu-shell 
scheme is considerably larger than the physical mass. This is somewhat unexpected. Let us 
recall that the naive treatment, which neglects the meson self energy, is the most used in 
practice. It has been known for a long time [p2| , p3| that the tu-shell scheme matter 
of principle, the correct procedure but on the basis of rough estimates it was assumed that 
neglecting the meson self energy would be a good approximation for the meson mass. We 
find here that this is not so. 

Regarding the consequences of removing the ghost, we find in table |I| that the effective 
parameters nis^v and gs,v in the no-ghost scheme are similar, within a few per thousand, to 
those of the naive scheme. This similarity reflects in turn on the predicted nuclear properties: 
the results shown in tables |T| and |^ for the no-ghost scheme coincide, within the indicated 
precision, with those of the naive scheme (not shown in the table). It is amazing that 
the outcoming parameters from such a sophisticated fltting procedure, namely the no-ghost 
scheme, resemble so much the parameters corresponding to the naive treatment. We believe 
this result to be rather remarkable for it justifles a posteriori the nowadays traditional 
calculations made with the naive scheme. 

The above observation is equivalent to the fact that the zero-momentum masses, m^,^, 
and the propagator-pole masses m^^^ are very similar in the no-ghost scheme. This implies 
that the effect of removing the ghosts cancels to a large extent with that introduced by 
the meson self energies. Note that separately the two effects are not small; as was noted 
above m„ is much larger than in the a;-shell scheme. To interpret this result, it will be 
convenient to recall the structure of the meson propagators. In the leading approxima- 
tion, there are three kinds of states that can be created on the vacuum by the meson flelds. 
Correspondingly, the spectral density functions p(g^) have support in three clearly separated 
regions, namely, at the ghost mass squared (in the Euclidean region), at the physical meson 
mass squared, and above the A^A^ pair production threshold (2M)^ (in the time-like region). 
The full meson propagator is obtained by convolution of the spectral density function with 
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the massless propagator (g^ + ie)~^ as follows from the Kalleen-Lehmann representation, 
eq. (0). The large cancelation found after removing the ghosts leads to the conclusion that, 
in the zero-momentum region, most of the correction induced by the fermion loop on the 
meson propagators, and thereby on the quadratic kernels K^p^), is spurious since it is due 
to unphysical ghost states rather than to virtual NN pairs. This can also be seen from 
figs. ID and ^ There, we represent the real and imaginary parts of Ks{p'^) respectively, in 
three cases, namely, before ghost elimination, after ghost elimination and the free inverse 
propagator. In all three cases the slope of the real part at zero momentum is equal to one 
and the no-ghost (sea 2nd) set of parameters from table | has been used. We note the strong 
resemblance of the free propagator and the ghost-free propagator below threshold. A similar 
result is obtained for the vector meson. 

One may wonder how these conclusions reflect on the sea energy. Given that we have 
found that most of the fermion loop is spurious in the meson self energy it seems necessary 
to revise the sea energy as well since it has the same origin. Technically, no such problem 
appears in our treatment. Indeed the ghost is found in the fermion loop attached to two 
meson external legs, i.e., terms quadratic in the fields. However, the sea energy used, namely, 
^sea _j_ ]^^^^ does not contain such terms. Quadratic terms would correspond to a mass term 
in E^'^ and a kinetic energy term in E^^, but they are absent from the sea energy due 
to the zero-momentum renormalization prescription used. On the other hand, terms with 
more than two gradients were found to be negligible [§. Nevertheless, there still exists the 
possibility of ghost-like contributions in vertex functions corresponding to three or more 
mesons, similar to the spurious contributions existing in the two-point function. In this case 
the total sea energy would have to be reconsidered. The physically acceptable dispersion 
relations for three or more fields have been much less studied in the literature hence no 
answer can be given to this possibility at present. 

V. SUMMARY AND CONCLUSIONS 

We summarize our points. In the present paper, we have studied the consequences of 
eliminating the vacuum instabilities which take place in the a-uj model. This has been done 
using Redmond's prescription which imposes the validity of the Kallen-Lehmann represen- 
tation for the two-point Green's functions. We have discussed possible interpretations to 
such method and have given plausibility arguments to regard Redmond's method as a non 
perturbative and non local modification of the starting Lagrangian. 

Numerically we have found that, contrary to the naive expectation, the effect of including 
fermionic loop corrections to the mesonic propagators (cj-shell scheme) is not small. However, 
it largely cancels with that of removing the unphysical Landau poles. A priori, this is a 
rather unexpected result which in fact seems to justify previous calculations carried out in 
the literature using a naive scheme. Actually, as compared to that scheme and after proper 
readjustment of the parameters to selected nuclear matter and finite nuclei properties, the 
numerical effect becomes rather moderate on nuclear observables. The two schemes, naive 
and no-ghost, are completely different beyond the zero four momentum region, however, and 
for instance predict different values for the vector meson mass. 

Therefore it seems that in this model most of the fermionic loop contribution to the meson 
self energy is spurious. The inclusion of the fermionic loop in the meson propagator can only 



16 



be regarded as an improvement if the Landau ghost problem is dcah with simultaneously. 
We have seen that the presence of Landau ghosts does not reflect on the sea energy but it 
is not known whether there are other spurious ghost-like contributions coming from three 
or higher point vertex functions induced by the fermionic loop. 



ACKNOWLEDGMENTS 



This work is supported in part by funds provided by the U.S. Department of Energy 
(D.O.E.) under cooperative research agreement #DF-FC02-94ER40818, Spanish DGICYT 
grant no. PB92-0927 and Junta de Andalucia grant no. FQM0225. One of us (J.C.) 
acknowledges the Spanish Ministerio de Educacion y Ciencia for a research grant. 



APPENDIX A: 

1. Meson self energies in the leading 1/N expansion 

As stated in the main text, the leading order in the 1/A^ expansion {N being the number 
of nucleon species) is achieved by considering 1-fermion loop and zero-boson loop Feynman 
graphs in the effective action. This corresponds to compute the meson self energies at 1-loop 
approximation. 

For the cr-meson, the bare self energy in terms of the Lagrangian coupling constant is 
obtained as 



^B,.(p';M,C,£) = -^e 



2e 



(27r)4-2^ 

i 



Tr 



J) + }/t- M + ie 



— M + 16 



Imposing zero-momentum renormalization we get 



47r2 [ V 

where the function Ia{y) is defined as 

Ia{y) ^ [ dx log[l -yx{l 
Jo 



- ^p' ] 4 



' p2 \ p2 



xj — ie\ 



'1- 






y 




- 1 







1-^+1 



-2 y <0 

- 2 < y < 4 

yrrilog^-^-2-zyrr|. A<y. 



(Al) 



(A2) 
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The cj-meson self energy is obtained in a similar way but taking care of its Lorentz 
structure, 

which is highly simplified by baryonic current conservation, 

^'bAp') = {-g'-" + ^) ^bAp') ■ (A4) 

The explicit expression of the cu-meson self energy renormalized at zero-momentum is 

n,(p2) = UbAp') - nB,.(o) - h'bMp' 



2. Poles and residues 



The relation between the Lagrangian mass of a boson, m, and its physical mass, mgh, is 
given by 



m, 



sh 



m 



Him, 



shJ 



0, 



(A6) 



where the self energy, n(p^), is assumed to be renormalized at zero momentum. From the 
expression of its self energy given above we find that the a-meson physical pole, m^, can 
be obtained in terms of the Lagrangian parameters Qg and by solving the transcendental 
equation 



m„ 



4 ml 

4 + 

3 M2 




1 arcsin 



rrin 



2M 



(A7) 



Similarly, the equation to solve for the uj particle is 



47r2 



A hml 2 ml 
3 9 M2 3 I M2 



4M2 m, 
— 1 arcsin 



m^ 



2M 



(A8) 



It is interesting to note that sometimes the combination = M^g'^/m'^ is taken to be fixed 
by nuclear matter properties. This allows one to write the Lagrangian coupling constant, 
g, as a function of C and the Lagrangian mass, m. Inserting the omega version of this 
expression into the previous equation permits to solve the Lagrangian mass in terms of 

are conveniently extended to the 



and the physical u mass, m^. If the eqs. (|A7|-p 
msh-complex plane they can be used to obtain the Landau ghost masses as well (better 
expressions for numerical calculation are found in the main text in eqs. (|35D and (|33|) ). 
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Once a Landau pole has been computed, the value of its zero-momentum residue, —Rg, 
is easily obtained as 

-Rg = 1- n'(-M2 ) . (A9) 

The particular expressions of this equation for the a and u meson are given in eqs. (0) and 
(P^) respectively. 
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TABLES 





9s 




rria 




rriy 




Tin 




Q Of,? 


44Q Q 


43Q 8 


1 3 81 


783 


673 6 




cij-shcll 




389 8 


37Q 8 


1 1 78 


Ql 8 


783 


sea 0th 


no-ghost 


5.996 


370.9 


368.3 


14.86 


978.6 


783 






(5.928) 


(368.8) 




(10.17) 


(786.1) 






naive 


5.922 


368.4 


365.9 


10.13 


783 


711.2 




a;-shell 


6.846 


425.9 


420.3 


11.78 


910.8 


783 


sea 2nd 


no- ghost 


6.664 


410.8 


406.5 


14.86 


978.6 


783 






(6.544) 


(407.1) 




(10.17) 


(786.1) 






naive 


6.536 


406.6 


402.6 


10.13 


783 


711.2 



TABLE I. Zero momentum renormalized Lagrangian parameters in several schemes. Masses 
are in MeV. The meaning of the labels no sea, sea 0th, sea 2nd, w-shell and no-ghost are as in 
table III. The naive scheme corresponds to not including the meson self energy. The numbers in 

brackets are the zero-momentum parameters of the effective action for the no-ghost scheme, rhg^v 
and gs,v In all cases, mg- and ra^ stand for the poles in the meson propagators after including the 
one fermion loop self energy and using the corresponding Lagrangian parameters. Note that by 
construction the vector meson parameters coincide in the sea 0th and sea 2nd cases. 





Rq 


Mg 


TDV 
JXq 


Mh 


sea 0th 
sea 2nd 


1.748 
1.584 


4605 
3863 


0.6090 
0.6090 


1457 
1457 



TABLE II. Residue (up to a sign) and mass (in McV) of the ghosts in the zero-momentum 
renormalized meson propagators using the no-ghost sets of Lagrangian parameters in table I. 
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B/A (MeV) 




no sea 


sea 0th 


sea 2nd 


exp. 






no-g 


host 


cj-sheU 


no-ghost 


cu-shell 




2oCa 


6.28 


6.00 




6.10 


6.33 


6.43 


8.55 


iNi 


7.24 


6.51 




6.60 


6.80 


6.90 


8.64 


907^ 


8.36 


7.99 




8.07 


8.22 


8.30 


8.71 




8.81 


8.43 




8.50 


8.62 


8.69 


8.36 


208pb 


9.84 


9.55 




9.61 


9.70 


9.76 


7.87 



TABLE III. Binding energy per nucleon of some closed-shell nuclei computed in several ways: 
not including the Dirac sea in the parameter fixing (no-sea), including the Dirac sea at lowest order 
in a gradient expansion (sea 0th), including the Dirac sea at all orders (sea 2nd) and the exper- 
imental values (exp.). The entry w-shell corresponds to use the set of parameters that reproduce 
the (j-meson mass after including the meson self energy. The entry no-ghost corresponds to the 
parameters obtained by applying Redmond's prescription. 











m.q.c.r. 


(fm) 








no sea 


sea 0th 


sea 2nd 


exp. 






no-g 


host 


a;-shell 


no- ghost 


a;-shell 




40r'„* 
2o*^a 


3.48 


3.48 




3.48 


3.48 


3.48 


3.48 


iNi 


3.72 


3.79 




3.79 


3.79 


6.80 




907^ 
40 


4.22 


4.23 




4.24 


4.25 


4.25 


4.27 




4.60 


4.66 




4.66 


4.68 


4.68 




208pb 


5.35 


5.39 




5.39 


5.41 


5.41 


5.50 



TABLE IV. Mean quadratic charge radii of several closed-shell nuclei. Meaning of the labels 
and experimental values as in table III. 
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FIGURES 



FIG. 1. Real part of the inverse scalar meson propagator Ksip^) as a function of the squared 
foTir momentum using the no-ghost (sea 2nd) set of parameters of table I. The dashed line represents 
the one-loop result without ghost subtraction. The solid line is the result after ghost elimination. 
The dotted line shows the free inverse propagator. In all cases the slope at zero momentum is 
unity. Units are in nucleon mass. 

FIG. 2. Imaginary part of the inverse scalar meson propagator KsijP')- Units and meaning of 
the lines as in fig. 1. 
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FIGURES 



-3 




FIG. 1. Real part of the inverse scalar meson propagator KgijP') as a function of the squared 
four momentum using the no-ghost (sea 2nd) set of parameters of table I. The dashed line represents 
the one-loop result without ghost subtraction. The solid line is the result after ghost elimination. 
The dotted line shows the free inverse propagator. In all cases the slope at zero momentum is 
unity. Units are in nucleon mass. 
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FIG. 2. Imaginary part of the inverse scalar meson propagator Ksip^)- Units and meaning of 
the lines as in fig. 1. 
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